Measuring microwave quantum states: tomogram and moments 
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Two measurable characteristics of microwave one-mode photon states are discussed: a rotated 
quadrature distribution (tomogram) and normally/antinormally ordered moments of photon cre- 
ation and annihilation operators. Extraction of these characteristics from amplified microwave 
signal is presented. Relations between the tomogram and the moments are found and can be used 
as a cross check of experiments. Formalism of the ordered moments is developed. The state purity 
and generalized uncertainty relations are considered in terms of moments. Unitary and non-unitary 
time evolution of moments is obtained in the form of a system of linear differential equations in con- 
trast to partial differential equations for quasidistributions. Time evolution is specified for the cases 
of a harmonic oscillator and a damped harmonic oscillator which describe noiseless and decoherence 
processes, respectively. 
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I. INTRODUCTION 

Non-classical states of light are of great interest not 
only from the viewpoint of academic science but also 
from the viewpoint of their promising quantum informa- 
tion applications such as cryptography, communication, 
etc. Any quantum experiment consists of the follow- 
ing procedures: preparation of a desired quantum state, 
time evolution of the state (its transformation through 
a quantum channel), and a measurement. The latter 
one is a keystone to get an insight into quantum world. 
The result of a measurement is usually a set of outcomes 
distributed in accordance with the quantum nature of 
the state. Given a numerous number of identically pre- 
pared quantum states and neglecting the memory effects 
of channels, one can perform many individual measure- 
ments and gather the statistical information. We refer to 
the obtained information as a measurable characteristic 
of the state. At this stage one encounters a problem how 
a measurable characteristic is related with the quantum 
state itself, i.e. how to represent the results of measure- 
ments. 

A quantum state is usually described by wave function 
or density operator Alternatively, different repre- 
sentations of the density operator are widely used, e.g., 
the Wigner function 0]. Recently, a probability repre- 
sentation of quantum states was introduced In the 
probability representation, quantum states are identified 
with fair probability distributions. The probability rep- 
resentation was introduced in connection with homodyne 
measuring the photon quantum states by means of opti- 
cal tomography. Using the homodyne detection scheme 
gives rise to the optical tomogram w(X, 0) which has a 
meaning of the probability distribution of a single ro- 
tated quadrature component Xg = q cos 6 + p sin in the 
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phase space of one light mode (see, e.g., the review Q). 
The aim of homodyne measurements [5| was to obtain 
the Wigner function identified with the quantum state. 
Then the reconstructed Wigner function was used to ex- 
tract an information on the physical properties. The 
measurements of this kind in optical frequency domain 
were fulfilled, e.g., in H,0|, and the detailed review of the 
experimental results is presented in Ref. . Though op- 
tical tomogram is a fair probability distribution, it was 
interpreted in all the experimental works as a techni- 
cal ingredient providing tool to obtain the Wigner func- 
tion. The latter one was interpreted as a 'real state'. 
On the contrary, following the ideas of the probability 
representation [3, H| (see also the recent review [l(| and 
the paper [HI), the optical tomogram and other kinds of 
tomographic-probability distributions like symplectic to- 
mogram [l2| were considered as a primary object contain- 
ing complete information on quantum state. In view of 
this fact, one does not need reconstructing any quasidis- 
tribution including the Wigner function. The reconstruc- 
tion procedure produces extra inaccuracies related to the 
useless elaborating the experimental data by means of 
Radon integral transforms. The experiment on direct 
checking purity-dependent uncertainty relations [HI, [l4| 
and measuring the photon state purity and temperature, 
without reconstruction of the Wigner function or another 
quasidistribution, was performed recently in Ref. [l5[ . In 
this experiment, the optical tomogram was measured and 
considered as a primary object determining the quantum 
state. Moreover, the above physical characteristics were 
expressed in terms of the tomogram, which is nothing 
else but directly measurable alternative to the density 
operator and the Wigner function. 

On the other hand, at optical frequencies there is an- 
other measurable characteristic of the state. Indeed, us- 
ing the heterodyne detection scheme gives rise to the 
Husimi function Q(q,p) which also contains the full in- 
formation about the quantum state Since both to- 
mogram and Husimi function are extracted from exper- 
imental data, the relation between them (see, e.g., jlfjj ]) 
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can be considered as a cross-check of the experiment ac- 
curacy. Closely related to the Husimi function are or- 
dered moments {(a*) n a m ) and (a k (a*) 1 ), where a 1 and a 
are photon creation and annihilation operators, respec- 
tively. In fact, the moments also contain the complete 
information on the state (see, e.g., [l7l - [2"ol | ) . An optical 
scheme for measuring moments is proposed, e.g., in |2lj . 
The moments are complex functions of integers. The 
measurements give Husimi function-type characteristics 
corresponding to a two-mode state, from which the mo- 
ments for one-mode state are calculated. In this sense, 
the method of moments is not completely direct mea- 
surement in contrast to the optical tomography. Also, 
the moments are not a probability distribution as the 
tomogram is. Nevertheless, if the moments are known, 
then all desired characteristics of the quantum state can 
be calculated. 

The importance of microwaves in up-to-date quan- 
tum technologies can scarcely be overestimated. Using 
microwave wavelengths transforms the scale of experi- 
ments as compared to standard optical ones. As a re- 
sult, at microwaves the highest-quality superconducting 
resonators [22| are built and successfully used in the mi- 
crowave version of cavity quantum electrodynamics. Be- 
ing applied in one-dimensional resonators together with 
superconducting qubits, microwaves have opened an op- 
portunity to achieve the strongest ever coupling between 
the electromagnetic field and an artificial atom (qubit) 
within compact and integrable electrical circuits [23j . Ex- 
perimentally realized couplings of microwaves to trans- 
mons [24| and mechanical oscillators [25} made them a 
really significant tool for further progress of quantum in- 
formation. 

On the other hand, microwave quantum states are of 
great interest per se as carriers of quantum information. 
However, measuring quantum state of electromagnetic 
field changes dramatically and becomes a challenge when 
microwave radiation is under investigation. Detection of 
microwave field (especially itinerant modes) is compli- 
cated by the low efficiency of photodetectors, although 
suggestions of high-efficiency microwave photon detectors 
are made [2(|[27j- As a result, it is unlikely to carry out 
the photon counting measurements. In order to register 
the radiation reliably, amplifiers are widely used [HI, [29| 
though they unavoidably add some noise. In addition, 
interferometry experiments in optics extensively utilize 
beam splitters, whereas in microwave engineering the role 
of a beam splitter for a single mode is played by hy- 
brid junctions or power dividers [30l [3l| . Similarly, a mi- 
crowave signal is mixed with a local oscillator microwave 
field (via a nonlinear device called mixer) instead of pass- 
ing a signal mode and a local oscillator mode through a 
beam splitter in optics. 

In spite of the challenges dealt with microwave ra- 
diation, the measurements of the 'optical' tomogram 
and ordered moments are feasible and were reported 
recently [32| - [35j . Really, using a homodyne detection 
scheme with phase sensitive amplifiers less noisy than 



high-electron mobility transistors (e.g., a Josephson para- 
metric amplifier |36| ) enables, in principle, measuring 
tomogram w(X, 9). This experimentally accessible dis- 
tribution w(X, 9) differs from the genuine tomogram of 
the state w(X, 9) because of the noise presented. How- 
ever, the information of the probed state can still be ex- 
tracted from the data. An analogue of the 8-port ho- 
modyne detection of optical photons is realized at mi- 
crowave level by phase-insensitive amplifiers and an in- 
phase quadrature (I/Q) mixer. The I/Q mixer provides 
two outputs q and p described by a single envelope func- 
tion S — q + ip. Statistics of experimentally measurable 
quantities q = ReS and p — ImS allows constructing the 
histogram Q(q,p). If there were no extra noise added, 
this quasi-probability distribution Q(q,p) would be an 
appropriately scaled Husimi function Q{q,p) of the quan- 
tum state. Averaging the Husimi function with complex 
function (S*) l S k = (q — ip) l (q + ip) k results in the mean 
value of anti- normally ordered operator (see, e.g., [13]), 
namely, (d fc (a^)'). 

The aim of our paper is to consider both methods 
of measuring microwave quantum states, viz., the ho- 
modyne tomography and the measurement of ordered 
moments. We briefly discuss how to extract moments 
from data corresponding to the amplified signal. Us- 
ing different experimental schemes, an access to nor- 
mally and antinormally ordered moments can be real- 
ized. The relations between the tomogram and the mo- 
ments are derived, which in case of 'optical' tomogram 
and normally ordered moments coincide with relations 
found in Refs. 38, 39]. This connections can be uti- 
lized to perform cross-check of measurements. We also 
suggest a test in the form of inequalities for experimen- 
tally measured moments (both quadrature and creation- 
annihilation ones). This test is equivalent to checking 
uncertainty relations but avoids unnecessary reconstruc- 
tion of the Wigner function or density operator. Finally, 
we cannot help considering time evolution of the quan- 
tum state in terms of measurable characteristics. Time 
evolution in the tomographic probability representation 
was considered in Refs. [3l. [Tollioj. In this paper, we fill 
this gap for moments. Namely, the Moyal equation [4l[ 
for the Wigner function is rewritten in terms of moments 
as well as an cigenstatc problem of Hamiltonian is for- 
mulated and non-unitary evolution of damped electro- 
magnetic field oscillator is considered for a partial case 
of generic study of the quantum oscillator with dissipa- 
tion (42|. The latter problem is instructive to clarify de- 
coherence phenomena in microwave experiments, where 
the decoherence occurs due to a finite conductivity of 
waveguide walls or a lossy dielectric. 

The paper is organized as follows. 

In Seaim normally and antinormally ordered moments 
are extracted from amplified signals and calculated for 
examples of Fock, coherent, even/odd coherent [H|, and 
thermal states. In Sec. IIIII a brief review of symplectic 
and optical tomograms is given. In Sec. IIV1 relations be- 
tween the tomogram and the moments are derived as 
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well as purity is expressed in terms of moments. In 
Sec. [V] inequalities for moments and generalized purity- 
dependent uncertainty relations are suggested in connec- 
tion with experiments such as that in Refs. [H, 13 EH > 
where a deterministic generation of microwave quantum 
states is demonstrated. In Sec. IVI| a unitary time evolu- 
tion and a problem of Hamiltonian eigenstates are formu- 
lated for moments. In Sec. IVII1 a particular case of the 
damped time evolution of moments is analyzed. Finally, 
in Sec. IVIIH conclusions and prospects are given. 



of linear equations of the form 



II. 



MOMENTS 



To anticipate, let us outline a linear parametric ampli- 
fier as a constituent of microwave engineering. The am- 
plifier changes the quantum state of the signal and idler 
incident modes by an SU(1,1) transformation (46|, i.e. 
the Bogoliubov transformation, which has the following 
form in the Heisenberg picture (see, e.g., (47l l4§j): 



b s = y/ga + \[g - lh\ 

h = \Jg - ia f + y/gh, 



(i) 

(2) 



where b s and hi are annihilation operators of the amplified 
signal and idler modes, respectively, a and h are annihila- 
tion operators of the original microwave signal and extra 
noise modes, respectively. Depending on the operation 
of the amplifier, the output mode can be described by 
either Eq. (TTJ) or Eq. ©. 

If the experiment exploits a heterodyne detection 
scheme (33j-[35|, the measured data represent nothing 
else but an appropriately scaled histogram of the Husimi 
function ff amp (q, p) , which differs from the Husimi func- 
tion of the original signal H(q,p). However, the informa- 
tion about the quantum state can be revealed by means 
of moments. In fact, averaging H am p(q,p) with com- 
plex function (S*) l S k = (q - ip) l {q + ip) k results in the 
mean value of anti- normally ordered operator b k (V*) 1 [4!| , 
where b is determined by either (TIJ or ([2]). For the sake 
of brevity, we will consider the case b = b s only, which is 
readily rewritten for the opposite case b = bi if needed. 

Thus, experimentally accessible moments read 



(*+j')/2 



(W) = EEu Mi i" 

i=0 j=o v / \ J 

X (g _ i)(*+'-*-j)/ 2 (a*(a t ) J ')<(ft t ) fc - < fe , - J '), (3) 

where the signal mode a and the noise mode h are as- 
sumed to be independent. Once gain g is known, average 
value of the normally-ordered noise operators ((fv) m h n ), 
m,n = 0,1,... can be measured by analyzing vacuum 
state |0) of mode a for which (0\a l (a^ ) J |0) = HSy. Sub- 
stituting this result in @, we obtain an infinite system 



(vac 



k\l\ 



min(k.l) 

§ i\(k-i)\(i-i)\ 



■xg\g- i) (fc+z)/2_i ((h t ) fc_i h i ~ i )- 



(4) 



Solution of this system is readily found for lower mo- 
ments. If n — or m — 0, we obtain 

(ft)°h n ) = (g- l)-"/ 2 (vac|S0(6t)«|vac), (5) 
(fork ) = (g- l)-'"/ 2 (vac|^(6t)°|vac). (6) 

Similarly, the cases n = 1 or m = 1 yield 

((tifk 1 ) = (s-ir^/VaclS^lvac) 
-ng(g - l)"^/^^®)"- 1 |vac), (7) 
((hXrh 1 ) = (g - l)-( ro+1 )/ 2 (vac|6™(St) 1 |vac) 
-mg(g - l)^(" l+1 )/ 2 (vac|6™- 1 (St)°|vac). (8) 

Finally, if moments up to the fourth are of interest, i.e. 
n + m < 4 as, e.g., in (35[, then we add the missing 
moment 



(vac|6 s 2 (6t) 2 |vac) 



(9 - 1) 

-4g( 5 -l)- 2 (vac|6 s 1 (St) 1 |vac) 
+2 ff 2 ( 5 -ir 2 (v aC |6°(6t)Vc). (9) 



Once normally ordered moments of noise mode are cal- 
culated, the antinormally ordered moments of the mi- 
crowave radiation mode are found from experimental 
data by solving the linear system ([3]) with respect to 
((^(aJp). In the same way, the choice b = b\ enables 
calculation of the antinormally ordered moments of noise 
mode and normally ordered moments of the original mi- 
crowave mode. Relations between differently ordered mo- 
ments are given in Ref. [Bol ]. 

Since ordered moments are extracted from experimen- 
tal data, it is reasonable to compare them with theoret- 
ical values. Let us consider examples of quantum one- 
mode states and calculate the corresponding normally 
and antinormally ordered moments. 

As far as Fock states \N) are concerned, we have 

N\ 

(N\(tf) n a m \N) = { (N-n)\ lin = m ^ N > ( 10 ) 

otherwise; 



(N\a k (tf) l \N) = ^±Ms k . 



kl- 



(11) 



Among mixed states, we point out a thermal state 
given by the unitless temperature T. Moments are obvi- 
ously related with the average photon number and read 



{(a ) a ) thermal = 



n\S n 



(a fe (at)<> 



thermal 



(eVT _ i) 
k\S kl 

'i - e-vn 



(12) 
(13) 
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Classical-like states are represented by a family of co- 
herent states | a), i.e. eigenstates of the annihilation op- 
erator a\a) = a\a). These states can be generated, e.g., 
by masers, and represent an electromagnetic field of a 
local oscillator if \a\ ^> 1. The ordered moments are 



(a\(a^) n a m \a) = (a*) n a T ' 



min(fc,Z) 



(a\a k (a)) l \a) = £ 



p=0 



fc \ I I 



P J \P 



p\a k ~ p (a*) l -P. (15) 



Other examples of nonclassical pure states of light are 
even/odd coherent states |±) = A/± (|a) ± | — a)) with 

Af± = [2(1 ± exp(-2|a| 2 ))] 1 , which were introduced 
in Ref. [43| . Transformation of these states after passing 
a linear amplifier is considered in Ref. [5l| . Here, we 
calculate all the ordered moments and then explore their 
non-unitary evolution in Sec. IVII1 We have 



(16) 



(±\(a?) n a m \±) =Nl{a*) n a m 

X [l + (-l)"+ m ± e" 2|Q|2 ((-1)™ + (-l) m ) 

min(fc,Z) 



{±\a\a^\±)=AflY. 



p=0 



P \P 



p\ a k ~ p (a*y- 



l + (-l) fc+z - 2p ±e- 2|Q|2 ((-l) fc - ;t, + (-l) i - p )| . (17) 



III. 



TOMOGRAM 



As it is stated in Sec.|H apart from the density operator 
p, the state of one-mode photon state can alternatively 
be described by the symplectic tomogram [l2j 



w e (X, p, v) = Tr [pS(X - pq- vp)] , 



where X is the real argument, p and v are real param- 
eters, q and p are quadrature operators. In what fol- 
lows, we will be interested in the measurable character- 
istic, viz., the so-called 'optical' tomogram w{X,6) — 
w s (X,p = cos9,v — sin#), where 9 € [0, 2ir] is the phase 
of the local oscillator. Here, we use the adjective 'opti- 
cal' as a conventional one, although microwave quantum 
states are under study. Unless otherwise stated, we will 
refer to the 'optical' tomogram as tomogram simply. 

Analysis of the microwave quantum states employs lin- 
ear amplifiers, so the tomogram w amv (X, 9) of the am- 
plified signal is only accessible. Our goal is to obtain 
the relation between the tomogram of amplified quantum 
state and the tomogram of the original one. As shown in 
Ref. [HI , the amplification ((T|) is represented by a convo- 
lution of the signal and idler fields in phase s pac e. The 
Wigner function of the amplified signal reads [52j 



1 



W amp (a) - J </"' i W ( ) n- 1Joi ,..( > ) ■ 



where a = (q+ip)/ \/2. Assuming g ^> 1, Eq. (fl9l) can be 
readily simplified. Integrating such a simplified Wigner 
function with delta- function S(X — q cos 9 — p sin ff) yields 
the convolution expression for tomograms 



1 



(14) WampPM) = —Jd 2 /3w\^ — 

xW noisc (/3). 



V2 



(20) 



Associating w amp (X, 9) with measured data, it is pos- 
sible, in general, to get insight into the state tomogram 
w(X, 9) itself by fulfilling the deconvolution of (|20[) . pro- 
vided the thorough knowledge of the noise state. The 
noise mode is usually assumed to be in thermal state 
with an experimentally controllable temperature T [34j . 
If this is the case, then 



^amp 



(X,9) 



x/2 



na 2 g 



w 



X 



Y, 9 exp 



2a 2 



dY, (21) 



where Gaussian variance a is expressed in terms of the 



unitless temperature T via formula u = y| coth ^ . 

While processing the experimental data, it is con- 
venient to deal with tomographic moments (Xg) — 
/_ ~ X r w(X, 9)d9. These moments are known to contain 
the full information about the quantum state. Using f|2 1 [) . 
we get the connection between tomographic moments of 
the amplified signal (Xg) amp and tomographic moments 
of the original microwave photon state (Xg): 



(18) (X r e ) amp = 



[r/2] 

r E 

1=0 



21 



(21 -mx; 



r-2l\ _2l 



(22) 



where (21 - 1)!! = 1 • 3 • 5 • . . . • (2/ - 1) and (-1)!! = 1. 
All the tomographic moments (Xg) up to a desired order 
R can be readily obtained from experimentally accessi- 
ble tomographic moments (Ag) am p, r = 0, 1, . . . , R, by 
inversion of formula (|2"2")l . After that, the tomographic 
probability-distribution function can be also expressed 
through tomographic moments. 



IV. RELATION BETWEEN TOMOGRAMS 
AND MOMENTS 



Let us now derive relations between measurable mo- 
ments and measurable tomographic distributions. These 
relations are to be used as a cross check of two approaches 
to probing microwave quantum states. 

Using the formalism of characteristic functions 



\a 



) and 



e ), it is not hard to prove 



(19) that the Wigner function is expressed through normally- 
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ordered and antinormally-ordered moments as follows: 



following inverse relation: 



W(a) = 



d 2 x x n {-x*y 



x exp 



-|A| 2 + aA* -a*\ 



W(a) = [d*\[J2 
\ k,i 



(23) 



(a k (tf) l )(-\*) k \ l 

7T 2 kU\ 



x exp 



|A| 2 + aA* - a* A 



(24) 



It is worth noting that the expression (|23|) implies in- 
tegration over complex plane and then summation on 
n, m = 0, 1, . . ., whereas formula ([Ml) avoids singularities 
if the summation is done before the integration. As it is 
shown later, this peculiarity leads to simpler expressions 
for normally ordered moments than for antinormally or- 
dered ones. In Ref. [53| a discussion is presented how 
to reveal negativity of the Wigner function (|23p by the 
analysis of higher-order moments. 

Our goal is the relation between the 'optical' tomogram 
and the moments but at first we exploit the known map- 
ping of the Wigner function onto symplectic tomogram 
(see, e.g., [II]) and obtain after simplification 



w(X, fi, v) — —= exp 



(p + iis) n (p- ivy 



x- 



/i 2 + V 1 

H 



E 



{(a)) n a m ) 



J.! 17?. 



n+m 



w(X, \i, v) 

E 



i 

2^ 



re? 



dt; exp 



k,l 



k+l 




(p 2 + v 1 ) - %ix 



where H^(X) is the Hermite polynomial of degree N. 
Substituting cos# for \i and sin# for v in (|2"5|) - (l2l)l) . the 
'optical' tomogram w(X, 9) is expressed in terms of or- 
dered moments 



w(X,9) 



V2 n + m n!m! 



2tt 



E 

k,l 



kin 



y/2 



(28) 



((a t ) n a m \ 



lll7?.l 



27rV2"+ m (n + to)! 

U 

v i(m—n)8 i i ( xr\ 



27T +00 

d0 y w(X, 6) 

— oo 

(29) 



27rv / 2 Il+m (n + to)! 



dfl e^ m -^ e {H n+m {X e )), (30) 



where (H n+m (Xg)) depends on 9 and is obtained from 
polynomial H n+m (X) by replacing X r —> (Xg) for all 
r = 0,1,..., n + m. Thus, formula (|3"0"|) establishes a 
relation between normally ordered moments and tomo- 
graphic moments. We must note that the relation (|29|) 
was found previously in the papers [HI, |39[ . 

As far as antinormally ordered moments are concerned, 
a direct use of Eq. (|28p is sophisticated so we resort to the 
connection between antinormally and normally ordered 
moments and exploit result of Eqs. (|21)1) - (|3"0"|) to obtain 



kill 



2tv +oo 



min(fc,Z) 



J d9 J dX w(X,9) e 



i(k-l)6 



-oo 



E 



km 



2 p Hk+i-2 P (X) 
^ p\(k + l-2p)l 

i(fe 



(31) 



p=0 



p!(fc + ; - 2p)! 



In order to perform cross check of the experimental 
data (tomogram vs. moments), one can also compare 
measured tomographic moments (Xg) with tomographic 
moments predicted by the measurement of normally or- 
dered moments ((d)) n d m ) : namely, 



(X r ) = £ __- V irr _ n _ mAl (^) n a m )e^- m)e . (33) 



n \ m \ ( r-n-m\\ 
r— n— m is even 



The latter formula is obtained by combining a definition 
of the tomographic moment (Xg), Eq. (|2"T|) . and the 
following integral 



+00 



A , ' J ff A r(X)e- A dA 



Since both tomogram and moments are measurable in 
microwave experiments, Eq. (|27l) is the basis for a cross- 
check of homodyne and heterodyne detection schemes. 

Let us now address the problem to invert formulas 
(|2"7|) - (|2"5|) . Using orthogonality of Hermite polynomials 
and orthogonality of trigonometric functions, we have the 



'irr 



27— iV ^ r-iV ^| 



if r — N is odd or N > r, 

(34) 



otherwise, 



which can also be written in the form of two-dimensional 
Hermite polynomial [54j . 
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A. Purity in terms of moments 

In this subsection, we attack the problem how to deter- 
mine purity of the state without unnecessary reconstruc- 
tion of the density operator but dealing with measured 
normally ordered moments. Just as the measurable to- 
mogram is used for fast and reliable calculation of the 
purity parameter (l5l | , here we derive the purity in terms 
of moments only. 

To start with, the overlap Trfpi^] between two density 
operators p\ and /3 2 equals to the overlap of the corre- 
sponding Wigner functions. In view of the relation 
it is readily seen that 



Tr[pip 2 ] = 



n.m.k.l 



(-l) m+k {n + k)\ ^ 

I II III "n+fe,m+Z 

n\m\K\l\ 
x((at)"a m ) 1 ((at) fe a') 2 . (35) 



For example, if one is interested in how close the pre- 
pared state is to the vacuum one, it is enough to calculate 

(0|/5|0)=E fc (-l) fc (fc!)- 1 ((« t ) fc « fc )- 

Finally, the state is thoroughly described by moments 
((a') n a m ) and its purity 7f is 

- (-l) m+k (n + fc)! /fzi\nzm\/rzi\kzi\ 
71 = 1 [Tm °n+k,m+l{{a ') a )(( a ') «)• 



n.m.k.l 



(36) 

Purity of the thermal state (|12l) is ^thermal = 
tanh(l/2T) so the effective temperature T c g of the elec- 
tromagnetic field can be probed directly via measured 
moments as T e g = 1/(2 arctanh-zr). 



V. UNCERTAINTY RELATIONS 

The role of uncertainty relations in quantum mechanics 
can scarcely be overestimated. Therefore, a direct prob- 
ing of uncertainty relations in experiments is of great in- 
terest. The information in amplified microwave signals is 
obscured by the noise. However, as it is shown in Sec-ITu 
the ordered moments of the original microwave mode can 
still be extracted from the data. The question arises itself 
whether such extracted moments satisfy the uncertainty 
relations. The violation of these inequalities would indi- 
cate the incorrectness of data processing or measurement 
of the data. 

Let us start with a conventional Schrodinger- 
Robertson inequality a qq a pp — a qp > 1/4. Written in 
the form of tomographic moments, i.e. ((AAg) 2 ) — 
(Xg) — (Xg) 2 , this inequality takes the form 55] 



((AX e f)((AX e+7T/2 f) 



((AA e+x/4 ) 2 ) 



-i(((AA e ) 2 ) + ((AA 0+x/2 ) 2 ) 



> i 07) 



and is to be satisfied for any angle 9 £ [0,27r]. Us- 
ing the established relation between tomographic mo- 
ments and normally ordered moments (|33j) . we have 
(X e ) = ((aV 9 + (a)e- 19 ) /y/2 and (X 2 ) = l/2+{tfa) + 
{{{a)) 2 )e i26 + (a 2 )e- i2e ) /2. Finally, we obtain the un- 
certainty relation in terms of moments 



((#&)-$)(&)) + a) -(tf)(&)) 2 

- «(at) 2 ) - (at) 2 ) ((a 2 ) (a) 2 ) > 0, (38) 

which turns out to be independent on 9. 

The inequality (|38| can be made stronger if we take 
into account the purity of the state. Indeed, we 
have [HI a qq a pp - a 2 qp > $ 2 (7r)/4, where <&{tt) « (4 + 
Vl6 + 97r 2 )/97r within the accuracy of 1%. Consequently, 
purity-dependent uncertainty relation in terms of nor- 
mally ordered moments is 

- (((at) 2 ) - (at) 2 ) ((a 2 ) - (a) 2 ) > (<f> 2 « - l) /4, (39) 

where the purity tt is to be calculated via moments in 
accordance with Eq. Q36p. Thus, we have formulated 
a self-consistent problem of checking purity-dependent 
uncertainty relations by using measurable moments only. 

Another drawback of inequality (|38| is that it exploits 
the lowest moments only, i.e. ((at)™a m ) with n + m < 2. 
Let us derive uncertainty relations in terms of moments 
such that they involve moments up to a desired order. In 
fact, suppose the operator F = z t+J2 n (yn^ n + z n(o^) n ) 
with arbitrary complex numbers zo,y n ,z n . Then the 
mean value (F'F) > for all zo,y n ,z n G C. This 
immediately implies that the quadratic form is positive- 
semidefinite. According to Sylvester's criterion, all lead- 
ing principal minors of the corresponding matrix are non- 
negative. For instance, we readily obtain the inequalities 
on moments up to the fourth order: 



/ (1) 



a) 



((at) 
V (a 2 ) 



(a) 


(at) 


(a 2 ) 


(fit a) 


<(« f ) 2 > 


(at a 2 ) 


(a 2 ) 


(act) 


(a 3 ) 


((at) 2 a) 
(a 3 ) 


((at)3) 
(a 2 at) 


((at) 2 a 2 
(a 4 ) 



((at) 2 ) \ 

((at) 3 ) 
(a(at) 2 ) 
((at) 4 ) 
(a 2 (at) 2 ) J 



> 0. 



(40) 



If normally ordered moments are measured, then one 
should replace (aat) = (a,ta) + 1, (a 2 a)) 
(a(at) 2 ) = ((at) 2 a) + 2(&), (a 2 ^) 2 ) 



ata 2 ) +2(a), 



((a 



4(a'a) + 2. Otherwise, i.e. if antinormally ordered mo- 
ments are experimentally determined, then one should re- 
place (ata) = (aat)-l, (ata 2 ) = (a 2 at)-2(a), ((at) 2 a) = 
(a(at) 2 )-2(at), ((at) 2 a 2 ) = (a 2 (fit) 2 ) -4(afit) +2. Note 
that formula (|38|> is nothing else but a condition on non- 
negativity of the second principal minor of the matrix 

gnj. 
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VI. UNITARY EVOLUTION AND 
EIGENSTATES 

In this section, we follow ideas of the seminal pa- 
Normally and antinormally ordered moments 



41] 



per 

are functions of non-commuting operators a and a< and 
can be considered as specific 'phase-space quasidistribu- 
tions'. Now we are aimed at deriving the laws which gov- 
ern the transformation with time of these 'phase-space 
quasidistributions'. In other words, the problem is to 
find time evolution equations for ordered moments. If 
this problem is solved, then such equations can be used 
as an alternative to the Schrodinger equation for the wave 
function, the von Neumann equation for the density op- 
erator, and the Moyal equation for the Wigner function. 
The crucial point is that ordered moments are experimen- 
tally measurable characteristics of microwave quantum 
states in contrast to wave functions, density operators, 
and Wigner functions. 

We start with correspondence rules between the oper- 
ators acting on the Wigner function and the operators 
acting on the ordered moments. Applying operators q, 
p, d/dq, and d/dp to the left-hand side of Eq. (|2U)) . we 
obtain 







(n) 



d_ 

dp 



mA 



-A 



(n) \ 
0,+lJ 

(n) > 
x 0, + lj 



(n) \ 
0,-1/' 

2V2V n£i -l,0 

to 



-mA 



(n) \ 

0,-lJ 



(41) 

(42) 

(43) 
(44) 



where a displacement operator A^ J for the normally or 
dered moments is introduced as follows: 



((tf) n+i a m+: >). (45) 



A^{(a^ n a m ) 



Arguing as above and employing Eq. ([24)) . the cor- 
respondence relations for antinormally ordered moments 
are found 



d 



d_ 

dp 

q O 
p 



J-/'A (a) + 

^2^+1,0 + 
i_(A (a) 



A (a) "I 



-A 



(a) • 
0,+l. 



IA 



(a) S 
0,-lJ 



■ZA 



(a) 



(46) 

(47) 

(48) 
(49) 



with a displacement operator A^ for the antinormally 



ordered moments being 

Ag(a fc (at)') := (d k +\d^+>) 



(50) 



The time evolution equation for ordered moments is 
readily obtained from the Moyal equation 



d 



d 



dt ^ P dq 



00 

^ (2/ 



(-iy d 2r+1 u d 2r+1 

l)!4 r dq 2 ^ 1 dp 2r + l 



by substituting normally or antinormally ordered mo- 
ments for W and using the correspondence table (|4ip- 
dB} or gSJ-dMl), respectively. 

The stationary Schrodinger equation, i.e. the eigen- 
state problem H\E) — E\E), transforms into the known 
equation for the Wigner function 



El 
2 



\_&_ 



(-l) r d 2r U d 



-)2r 



{2r)W dq 2r dp 



2r 



W E = EW E (52) 



and then to the equation for ordered moments according 
to associations and (|4T?)) -(l4T)l). 

In order to demonstrate the equations for moments, 
we consider the simplest case of the free evolution of the 
electromagnetic field, which is effectively governed by the 
harmonic oscillator potential U = q 2 /2. If this is the 
case, the lowest-order derivatives are only presented in 
Eqs. (|5ip and (I52[) . For normally ordered moments we 
have 



d_ 

dl 



i(n — m) 



A 



to 

+ + 1 



2 



((at)«a m ) = 0, (53) 
{{a)Td m ) E = E{{tf) n a m ) E . (54) 



It is not hard to see that the moments (TTU|) multiplied 
by e l(n ~ m)t are solutions of both Eqs. ([53]) and §4$ if 
E = N+ 1/2. 

Likewise, the antinormally ordered moments (|11[) mul- 
tiplied by e l ^ l ~ k ^ t are solutions of equations 



| + <(*-*) 



(a fc (&t)') = 0, 



(55) 



A (a) _ k+l + 1 
^+1, + 1 2 



(a k (a}) l ) E = E(a k (tf) l ) E . (56) 



VII. DAMPED EVOLUTION 

In practice, microwave transmission lines have loss due 
to a finite conductivity of waveguide walls or lossy dielec- 
tric (3l| . In many experiments the loss may be neglected. 
On the other hand, quantum superpositions are very vul- 
nerable to the relaxation and decoherence while interact- 
ing with the environment. In view of this, the effect of 
loss is of interest. 

To describe a microwave quantum state in lossy envi- 
ronment we make use of a damped quantum oscillator 
model [42|, [56[. The Wigner function obeys the linear 
Fokker-Planck equation of the form [56| 



d_ 

dt 



d 



p- 



' d 
^d~p p + 



d ' 

^ dp 


W 






M 


d 2 


d 2 \ 


7 2 d 2 ' 


2w V 


dq 2 ' 


v dp 2 ) 


u> dpdq 



W, (57) 



W — (51) where 7 is the damping coefficient and 10 = yl — 7 2 . 

In view of correspondence relations (l4"Tj) - (|4"4"l) . Eq. ([5T)l 
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transforms into the following equation for measurable 
normally ordered moments: 



9 -I > 
— — n — m 

dt v ; 



((at)"a™) 



n + m 



-Kl + iw-^Mn-ljALto 

-|(1 - icj- 1 7)m(m - l)A^ n l 2 l ((at)"a m ). (58) 



In the same way, one can construct a time evolution equa- 
tion for antinormally ordered moments. 

The infinite system of equations (1581) is linear. Thus, 
the partial differential equation (|5T[) is reduced to a lin- 
ear system of the first order differential equations. This 
means that a formal solution of the system is the matrix 
exponent which is easy to compute if one is interested 
in time evolution of the lower-order moments. More- 
over, any linear dynamics of the system (not necessar- 
ily quadratic as in Eq. ([571) ) can be predicted numeri- 
cally in terms of moments, whereas the dynamics of the 
Wigner function would require solving finite-difference 
equations with artificial mesh spacing. Conversely, mo- 
ments ((a') n a m } are endowed by a natural 'mesh spacing' 
An, Am = ±1 and, last but not least, are experimentally 
accessible. 

An example of non-unitary evolution of the normally 
ordered moments is presented in Fig.[TJ The odd coherent 
state |— ) (|16p with a = 0.5 evolves in time according to 



the system of equations where damping coefficient 
7 = 0.1. Observation of the moments at successive times 
provides snapshots of the decoherence process (cf. [57|). 

VIII. CONCLUSIONS 

To resume we point out the main results of our paper. 

We have considered two approaches to measuring mi- 
crowave quantum states, namely, the homodyne detec- 
tion scheme with the 'optical' tomogram as output and 
the heterodyne detection scheme with output in the form 
of ordered moments of photon creation and annihilation 
operators. A microwave one-mode quantum state can be 
identified cither with the tomogram w(X, 9); or the to- 
mographic moments (Xg), r = 0, 1, . . .; or the normally 
ordered moments ((a^)™a m ), n, m = 0, 1, . . .; or the anti- 
normally ordered moments (a fe (a^) z ), k,l = 0,1,.... It 
was shown how to extract these quantities from mea- 
surable characteristics of the amplified microwave signal. 
We suggest using the established relations between the 
tomogram and the ordered moments (|27|) — (f29|) . (|3Tj) as 
well as the relations between the tomographic moments 
and the ordered moments (|30l) . (pJ5| as a cross check 
of the experimental results obtained in Ref. [32[ and in 
Refs. [Mli|. 

As the normally/antinormally ordered moments are 
measurable and determine a quantum state, an effort to 
obtain new results for the ordered moments has been 
made. Indeed, purity is expressed through moments 
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and is used in purity-dependent uncertainty relation in 
terms of moments. Another result is that the moments 
are to satisfy a generalization of the inequality ([40)) . 
We have obtained the time evolution equation in terms 
of moments, which is informationally equivalent to the 
von Neumann equation for the density operator and the 
Moyal equation for the Wigner function. The energy level 
problem and the non-unitary evolution of the damped 
microwave electromagnetic field are also considered in 
terms of moments. The damped evolution is described 
by a system of linear differential equations on moments, 
which is beneficial (from the viewpoint of numerical anal- 
ysis) in comparison with a partial differential equation on 
the Wigner function. 

Since normally /antinormally ordered moments are 
of great interest, a construction of the star-product 



scheme [58] for such moments is a problem for further 
investigation and will be considered elsewhere. 
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